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ABSTRACT 

We present near-infrared imaging spectroscopy of the strongly- lensed z — 2.00 galaxy SDSS 
J120601.69+514227.8 ('the Clone arc'). Using OSIRIS on the Keck 2 telescope with laser guide 
star adaptive optics, we achieve resolved spectroscopy with 0.20 arcsecond FWHM resolution in the 
diagnostic emission lines [Om], Ha, and [Nil]. The lensing magnification allows us to map the velocity 
and star formation from Ha emission at a physical resolution of ~ 300 pc in the galaxy source plane. 
With an integrated star formation rate of ~ 50Mq yr _1 , the galaxy is typical of sources similarly 
studied at this epoch. It is dispersion-dominated with a velocity gradient of ~ ±80 km s" 1 and aver- 
age dispersion a = 85 km s _1 ; the dynamical mass is 2.4 x 10 Mq within a half-light radius of 2.9 
kpc. Robust detection of [Nil] emission across the entire OSIRIS field of view enables us to trace the 
gas-phase metallicity distribution with 500 pc resolution. We find a strong radial gradient in both the 
[Nn]/Ha and [Oni]/Ha ratios indicating a metallicity gradient of —0.27 ±0.05 dex kpc^ 1 with central 
metallicity close to solar. We demonstrate that the gradient is seen independently in two multiple 
images. While the physical gradient is considerably steeper than that observed in local galaxies, in 
terms of the effective radius at that epoch, the gradient is similar. This suggests that subsequent 
growth occurs in an inside-out manner with the inner metallicity gradient diminished over time due 
to radial mixing and enrichment from star formation. 

Subject headings: galaxies: high-redshift — galaxies: evolution — gravitational lensing: strong 



1. INTRODUCTION 

Considerable progress is being made with integral field 
unit (IFU) spectrographs on large telescopes in unravel- 
ing the internal properties of star-forming galaxies over 
1 < z < 5. During this period, the mass assembly of 
galaxies proceeds at a rapid pace and such measures 
promise valuable insight into the physical mechanisms by 
which young galaxies grow. Ear ly effort has focused on 
the dynamical properties (e.g. |Forster Schreiber et al. 
2009|)Law et alj|2009t f Jones et al.l 120101: iGnerucci et al. 
2010D . These studies have revealed systems which are 
dispersion-dominated with varying degrees of ordered ro- 
tation consistent with sources which may develop stable 
disks and central stellar bulges. The turbulent motions 
may be associated with vigorous star-formation driven 
by cold molecular gas accrete d along filaments fro m the 
nearby intergalactic medium (|Tacconi et al.l l2010). The 
highest resolution data, secured by coupling the angular 
magnification of strong gravitational lensing with laser- 
assisted guide star adaptive optics (LGSAO), has re- 
vealed the sizes of star-forming regions in z ~2-3 galax- 
ies suggesting that star formation is primarily triggered 
by gravitational insta bility rather than external merg- 
ers (Stark et al.l 120081: Pones et aLlfeOltt ISwinbank et al.l 
12010) . 

Attention is now shifting to understanding the evo- 



lution of the mass-metallicity relation (jTremonti et al.l 
2004). Metallicity is a key parameter which gauges the 
baryonic material already converted into stars thus of- 
fering insight into feedback processes proposed to regu- 
late star formation. The trend is usually explained via 
star-formation driven outflows, e.g. from energetic su- 
pernovae, which have a larger effect in low mass galax- 
ies with weaker gravitational potentials. However, other 
effects may be responsible at high redshifts where star 
formation timescales and feedback processes and their 
mass dependences differ. Much effort has been invested 
in measuring evolution in the relation using the inte- 
grated light fro m sa mples beyond z ^ 2 (lErb et al.ll2006l; 
Hallid av et alJ200a IMannucci et~al1l2009HHavashi et all 
2009). These pioneering surveys have shown that metal- 
licity decreases with redshift for a fixed stellar mass. 

A logical next step in developing a picture of the as- 
sembly history of disk galaxies would be to extend the 
kinematic progress made with IFU spectrographs to re- 
solved measures of the gas phase metallicity. This would 
enable the tracking of metallicity gradients as a func- 
tion of redshift and provide a direct test of models of 
star formation and assembly induced by cold gas inflows 
(|Brooks et al.ll2007D . However, even with LGSAO, the 
spatial sampling for representative z ~2-3 galaxies is >1 
kpc and insufficient for this purpose. Prior to large tele- 
scopes with improved diffraction-limited performance, 
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gravitatio nally-lensed sourc es offer a natural starting 
point (c.f. iStark et al.ll2008[ ). In this Letter we examine 
the practicality of this approach via a detailed study of 
the gravitationally-lensed arc SPSS J120602.0 9+51422, 
referred to as the 'Clone' arc (|Lin et al.ll2009t ). Hubble 
Space Telescope (HST ) images and spectroscopic follow- 
up bv lLin et al.l ()2009() confirm this to be a z=2.001 star- 
forming galaxy with intrinsic half-light radius 77, = 2.9 
kpc and SFR ~ 50M Q yr _1 , magnified an areal factor 
~30 by a foreground z=0.422 group of galaxies. Further 
considera tions of the len s ing co nfiguration have been pre- 
sented by Vcg etti et al.l (|2010t ). The object was selected 
for detailed study after considering an integrated near- 
infrared emission line spectrum using NIRSPEC on Keck 
2 (Hainline et al 2009). 

Throughout this paper we adopt a A CDM cosmology 
with H = 70 km s- 1 Mpc" 1 , ft M =0.30 and ft A =0.70. 
At z=2.00, 0.1 arcsec corresponds to 840 pc and the age 
of the universe was 3.2 Gyr. All magnitudes are in the 
AB system. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of SDSS J120601.69+514227.8 (the 
Clone) were taken with the O H-Suppressing Infra - Red 
Imaging Spectrograph (OSIRIS; lLarkin et al.l 120061 ) and 
the LGSAO system on the Keck 2 telescope on 19 May 
2010. The weather was clear with strong winds and the 
seeing was 1.5 arcsecond FWHM. Nonetheless, tip/tilt 
correction yielded a respectable Strehl of 0.15 at 2 /jm. 
We took exposures of 900 seconds each, dithering by 1.8 
arcseconds in order to achieve good sky subtraction while 
keeping the arc within the OSIRIS field of view. We ob- 
tained 4 exposures in the Knl band (covering Ha and 
[Nil]), 2 exposures in Hnl ([Oin]) and 3 exposures in 
Zn4 ([On]). The physical location of the OSIRIS field is 
shown on a HST color composite in Figure [T] 

Dark and bias subtraction, spectral extraction, wave- 
length calibration and construction of the relevant data 
cubes was accomplish ed using the OSIR IS data reduc- 
tion pipeline (ODRP; lLarkin et al.l l2006). Sky s ubtrac- 
tion w as done with the IDL code described in iDavied 
(|2007|) using adjacent pairs of images as sky reference 
frames. The final data cubes were stacked using a clipped 
mean and show strong [Oin], Ha, and [Nil] emission. 
The brightest individual pixels have a signal-to-noise of 
~10 in Ha, 5 in [Oin], and 2 in [Nil]. The [On] doublet 
was not reliably detected and H(3 lay in a region of low 
atmospheric transmission. 

The infrared standard star FS 133 was observed with 
the same instru ment configuration. Prec ise JHK UKIRT 
photometry from lHawarden et al.1 (120011) was m atched to 
the G8v stellar spectrum from IPicklesI (1998) to deter- 
mine the conversion from ADU to flux as a function of 
wavelength. We tested this calibration by comparing the 
ratio of observed flux (in ADU s _1 ) to the reference spec- 
trum (in ergs -1 cm -2 ) with an atmospheric transmissio n 
model generated by the ATRAN program (ILord II1992T) . 
The data reproduced all major atmospheric absorption 
features to better than 10 % accuracy. We checked fur- 
ther for systematic errors by comparing the flux ratio 
of bright emission lines in the integrated OSIRIS spec- 
trum of image A3 (Figure [J) with in dependent results 
from NIRSPEC (jHainline et al.ll2009f) . The results were 
consistent to within the uncertainties (6%). Overall we 



deduce a systematic uncertainty in flux ratios of ~ 10 % 
for lines observed with different instrumental configura- 
tions (such as [Oni]/Ha), and negligible systematic error 
for lines observed simultaneously ([Nil] /Ha). 

3. GRAVITATIONAL LENS MODEL 

Our OSIRIS field samples several multiple images of 
the z=2.00 galaxy (A3-5; Figure [1]) enabling us, in prin- 
ciple, to derive the source plane properties independently 
via each image. However, accurate modeling of the fore- 
ground gravitational lens is necessary to re alize this ad- 
vantage. We use the l enstool program ()Kneib et al.l 
119931 : Uullo et al.l [2007) to parameterize the mass dis- 
tribution of the lens, using the image positions as con- 
straints. 

Our lens model comprises 1 group-scale dark mat- 
ter halo, 3 galaxy-scale dark matter halos associated 
with G2, G3, and G4, and an external shear compo- 
nent. Galaxies Gl, G2, and G3 are parameterized with a 
singular isothermal ellipsoid (SIE). G4 is parameterized 
with a dua l Pseudo-Isothermal el liptical mass distribu- 
tion (dPIE; lEhasd6ttir et al.ll2007D with an assumed tidal 
radius of 0.68" (jVegetti et aLlkoiOD . The fit is very good, 
with a mean residual error of 0.2" in the image plane and 
0.06" in the source plane. The recovered parameter val- 
ues for this model are reproduced in Table 1. We find 
typical linear magnification factors of 1.05 x 12.5 (orthog- 
onal and parallel to the arc) for image A3 and 1.01 x 6.7 
for A4 with an uncertainty of 5%. The area l magnifi- 
cation , 28.1 ± 1.4, is in good agreement with iLin et al.l 
(120091) . 

Lenstool defines the necessary transformations to the 
source (intrinsic) plane for each image, giving the source 
plane position for each spatial pixel of OSIRIS. We use 
the IDL functions triangulate and trigrid to inter- 
polate the irregularly-sampled source plane data onto a 
uniform grid. The source plane resolution of 0.15 x 1.6 
kpc FWHM was measured by reconstructing the tip/tilt 
reference star in the source plane at various locations. 
Since we will smooth the data to improve the signal/noise 
we note here that the resolution increases to 0.3 x 2.3 kpc 
and 0.5 x 3.0 kpc when the image plane data is convolved 
with a Gaussian of FWHM = 0.25 and 0.40 arcseconds, 
respectively. 

4. KINEMATIC PROPERTIES 

Considering the kinematic data for Ha in the image 
plane where the instrument resolution and sensitivity 
are nearly uniform and easily characterized, the inter- 
nal properties of the Clone are fairly typical of those 
z ~2-3 galaxies for which data of comparable reso lution 
is available (e.g. IStark et al.ll2008HJones et al.ll2010D . Af- 
ter smoothing with a Gaussian kernel of FWHM = 0.25 
arcsec we determine the Ha emission line velocity and 
width from a Gaussian fit to the line profile. The intrin- 
sic velocity dispersion is determined by subtracting the 
instrumental resolution of <Ji ns t = 52 km s _1 (measured 
from OH sky lines) in quadrature from the best-fit line 
width. Only fits with signal-to-noise > 5 in Ha are kept. 

In Figure [2] we present source-plane velocity and dis- 
persion maps constructed from the A3 image (see Fig- 
ure [TJ. The 2D map and ID velocity profile both show 
a large gradient of Av = ±80 ± 18 km s _1 across the 
galaxy with some degree of asymmetry. The recon- 
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FIG. 1: (Left) Hubble Space Telescope color composite image of the Clone arc. Al-5 represent multiple images of the z=2.00 
source; the critical curve is shown in red. The OSIRIS pointings are indicated by the two green rectangles offset by 1.8 arcseconds 
in the East-West direction. Foreground lensing galaxies are labeled as Gl-4. (Right) Distribution of the key emission line fluxes 
in multiple images A3-5 in units of 10~ 18 erg/s/cm 2 (see text for discussion of optimum sampling). North is up and East is to 
the left. The critical line passes through the A3 and A4 components such that a small part of the arc (including A5) is imaged 
4 times in the OSIRIS field. The remainder of A3/A4 is imaged twice in the OSIRIS field. 
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TABLE 1: Best fit parameter values and la uncertainties for each component of our lens model. Components G2, G3 and 
G4 are assumed to be spherical mass distributions located at the galaxy centroids. The central position of Gl is in arcseconds 
relative to the galaxy centroid. 



struction of images A4-5 gives a consistent result with 
Av = ±86 ± 19 km s -1 . The velocity dispersion is typ- 
ically a = 130 ± 20 km/s in the center and shows sim- 
ilar asymmetry. The mean dispersion of individual pix- 
els is a — 85 ± 2 km s _1 , consistent with a = 80 ± 4 
km s _1 derived from the Hainline et al.l ([2009) inte- 
grated spectrum. This dispersion gives a dynamical mass 
Mdun = 5r h a 2 /G = 2.4 x 10 10 M Q for r h = 2.9 kpc 
(Hai nline et al.ll2009ft . 

The above observations are characteristic features of 
star-forming galaxie s observed at z ~ 2 (e.g. the SINS 
sample discussed bv lForster Schreiber et alT l2"009). The 
SINS galaxies have a median dispersion of 0{ nt = 130 
km s _1 and Av/<Ti nt — 0.4 — 2.0 for those classified as 
likely disk progenitors. Even those SINS galaxies with 
the most prominent rotational motion show asymmetric 
features, with model re siduals of > 100 km/s in both 
velocity and dispersion (jCresci et al.ll2009l ). 

5. GAS PHASE METALLICITY MEASUREMENTS 

Reliable metallicity measurements require knowledge 
of the electron temperature, ionization state, and red- 
dening, for which several emission lines must be observed. 
In practice, all metallicity measurements at high redshift 
rely on strong emission lines which can be detected in 
a reasonable integration time. Ratios of strong emis- 
sion lines are then used to infer the gas-phase oxygen 



abundance, calibrated locally using either direct mea- 
surements of the electron temperature or via photoioniza- 
tion models. One of the more po pular diagnostic s is the 
[Nll]A6584/Ha ratio (N2) which iPettini fc Pagell ()2004f l 
showed is strongly correlated with the oxygen abundance 
via: 

12 + log O/H = 8.90 + 0.57 x log [Nil] /Ha 

with an intrinsic la scatter of 0.18 dex. 

The main drawbacks of the N2 method are that 
[Nil] /Ha does not necessarily trace the global metallicity 
when (i) AGN and shock excitation contribute to [Nil] 
emission; (ii) secondary production of nitrogen leads to 
variation in the N/O ratio; and (hi) [Nil] cooling satu- 
rates (at 12 + log O/H > 9.0). The advantage of N2 lies 
in the proximity of the two required emission lines, such 
that systematic errors due to reddening and instrumental 
effects are negligible. 

Fortunately, we can use [Oiii]A5007 as an additional 
constraint on shock ionization, AGN, and variations in 
the N/O ratio. We smoothed the data cube spatially 
with a 0.4 arcsec FWHM Gaussian, providing a signal-to- 
noise > 3 for [Nil] across the entire OSIRIS field of view. 
We then fit a Gaussian profile to the Ha line and deter- 
mined the best-fit [Nil] and [OIII] fluxes for the same 
profile, requiring the Ha signal-to-noise to be > 10 (Fig- 
ure[TJ . We found no significant offset in the velocity when 
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FIG. 2: Kinematics of the Clone. Left: velocity and velocity dispersion maps in the image plane (top) and source plane 
(bottom); North is up and East is to the left. The critical line and corresponding caustic are overlaid on the velocity maps. 
Right: one-dimensional velocity (circles) and dispersion (triangles) profiles in the source plane extracted along a slit passing 
through the centroid of Ha emission and oriented to maximize the velocity gradient. The same velocity gradient is seen in both 
images (A3-4) of the galaxy, with redshifted gas at the center of the IFU and blueshifted gas at both edges of the field. 




FIG. 3: Source plane reconstruction of the Ha signal-to- 
noise, [Nil] /Ha, and [Om]/Ha maps from the A3 region of the 
arc. Axes are in kpc. The maps show a maximum [Nil] /Ha 
and minimum [Om]/Ha near the center of the galaxy, slightly 
offset from the Ha emission peak. 

fitting these lines independently. For each pixel we infer 
the [Oin] A5007/H/3 ratio (assuming Case B recombina- 
tion) and measure the [Nil] /Ha ratio ad opting an aver- 
age s elective extinction of E(B-V) = 0.28 (|Hainline et al.l 
2009). The resulting data show no evidence of shock 
ionization or AGN in the standard BPT diagnostic di- 



agram (Bal dwin et al.lfTMl . We can examine possible 
variations in the N/O ab undance using the photoi oniza- 
tion models presente d in IKewlev fc D opita (2002). The 
[Oii]/[Oiii] flux from lHainline et all (|2009D . corrected for 
reddening via the Ha/H7 ratio, suggests an ionization 
parameter \og(U) — —2.9 for the integrated spectrum. 
We compute the metallicity of each pixel from both the 
[Nil] /Ha and [Nii]/[Oiii] ratios using this ionization pa- 
rameter and find excellent agreement, with an average 
offset 0.02 and la scatter of 0.08 in 12 + log O/H. 

The source-plane reconstruction is done independently 
for the A3 and A4/A5 regions, yielding two indepen- 
dent measurements of the metallicity distribution. The 
results are shown in Figures |3] and 0J [Nil] /Ha peaks 
in the center of the galaxy roughly coincident with the 
global minimum of [Oni]/Ha, with [Nil] /Ha decreas- 
ing and [Oin]/Ha increasing at larger radii. This is 
the expected signature of a metallicity gradient as seen 
in local disk ga l axies (e.g. iMagrini et al.l [20071 ) . Using 
iPettini fc Pagell (|2004h 's calibration to convert the mea- 
sured [Nil] /Ha to an inferred metallicity 12 + \ogO/H, 
we extract pixels within a slit oriented along the direc- 
tion of highest magnification to maximize spatial sam- 
pling and resolution. We bin pixels within the extrac- 
tion slit based on their radius from the center of the 
galaxy, determined from the centroid of the global min- 
imum [Oin]/Ha. A linear fit to each set of data gives 

dlos d ° i /H = -0.31 ± 0.07 and -0.27 ± 0.08 dex/kpc, 
with central values of 12 + \ogO/H = 8.59 ± 0.04 and 
8.49±0.04 respectively. Taken together, we infer a metal- 
licity gradient of _ _o.27±0.05 dex/kpc. Vary- 
ing the central position by up to 800 parsec results in 
gradients typically shallower by 0.05 ± 0.05 dex/kpc, but 
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FIG. 4: Intrinsic metallicity gradient of the Clone arc de- 
rived using the N2 index. Measured values of [Nil] /Ha are 
shown on the right axis with inferred oxygen abundance on 
the left. Filled circles and crosses correspond to the source 
plane properties of images A4 and A3, respectively. The 
dashed lines are linear fits to the two data sets. Thick er- 
ror bars represent only the statistical uncertainty in [Nil] /Ha 
for each bin, while thin error bars include an additional un- 
certainty from the N2 index calibration as discussed in the 
text. 

a consistent central metallicity. The metallicity inferred 
for the global spectrum is 12 + log O/H = 8.42 ± 0.18 
(adop ting the scatt e r in t he N2 calibration), consistent 
with Ha inline et al.l (|2009ft within the uncertainties. 

The preceding analysis assumes that scatter in the N2 
calibration applies to individual measurements but not 
their differences. We also measured the metallicity gra- 
dient with scatter in the N2 calibration treated as an 
additional uncertainty in each measurement (and their 
differences), but with no overall systematic offset. The 
resulting central metallicity is 12 + logO/iJ = 8.53±0.06 

and the gradient is dl ° s d % /B = -0.27 ± 0.11 dex/kpe. 
As discussed above, the [Om] measurements suggest no 
significant error in the inferred metallicity from varying 
ionization parameter and N/O ratio, implying that scat- 
ter in the N2 calibration does not affect relative abun- 
dances such as g radients. This is supported by studies o f 
nearby galaxies (jRupke et al.ll2010t iKewlev et al.ll2010D . 
hence we adopt an uncertainty of ±0.05 dex/kpe in the 
metallicity gradient. Regardless, Figure |4] clearly shows 
a significant radial gradient with correlation coefficients 
of —0.98 and —0.89 for images A4 and A3 of the arc, 
respectively. 

6. DISCUSSION 

The Clone arc is evidently representative of UV-bright 
star-forming galaxies at z ~ 2 both in its integrated star 
formation rate of ~ 50A/q yr _1 and its resolved kinemat- 
ics. Via the detailed study of 3 emission lines, we have 
found a strong radial gradient in both the [Nil] /Ha and 
[Om]/Ha flux ratios, consistent with a metallicity gra- 
dient which we measure to be dX ° & d ^ H = —0.27 ± 0.05 
dex/kpe. This is stronger than that in local disk galax- 
ies, which typically h a ve gradients o f ~ — 0.05 dex/kpe 
(jBresolin et all 120101 : iMoran et all 120101 : iRupke et all 
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[2010HCarrera et al.ll200l iMagrini et al.ll2007[ h although 
the total range A log O/H ~ 0.5 is similar. The near- 
solar central metallicity 12 + log O/H = 8.5 ± 0.2 is 
comparable though somewhat lower than in local mas- 
sive disk galaxies. We can place a rough upper limit 
on the age of the Clone by considering the timescale 
t sf = i^ph < gpn = 500 Myr, suggesting that strong 
central enhancement occurs very rapidly. 

Negative metallicity gradients suggest that the inner 
regions are older and more gas-depleted than the out- 
skirts. The extremely steep gradient and similar range 
of metallicity in such a young galaxy as the Clone sug- 
gests an inside-out growth scenario in which the cen- 
ter is enriched rapidly while stellar mass and metallicity 
gradually increase at large radii. This is further sup- 
ported by kinematic evidence: the large velocity disper- 
sion should cause the metallicity gradient in the Clone 
to flatten due to radial mixing on timescales of order 
Tdyn = d/o = 70 Myr. Infall of metal-poor gas onto the 
outskirts of the galaxy is therefore necessary to maintain 
a strong gas-phase metallicity gradient, which will natu- 
rally flatten toward values observed in local disks when 
accretion stops. 

The chemical and dynamical properties of the Clone 
are consistent with cosmological simulations which show 
that accretion of cold gas onto a t urbulent disk leads to 
a bulge plus thin d isk morphology (|Ceverino et al.|[20Tol : 
iBrooks et al.l [2009). The "clumpy disk" phase can last 
for several 10 8 years, during which radial abundance gra- 
dients would be affected by the competing processes of 
radial mixing (which flattens the gradient) and infall of 
metal-poor gas (which steepens the gradient), ultimately 
keeping the abundance spread relatively similar. Mean- 
while the galaxy grows in size and stellar mass, while 
centrally-migrating clumps form a bulge which stabilizes 
the disk. 

As emphasized in Uones et al.1 (|2010l ) our results 
demonstrate the considerable advantage of using lensed 
galaxies to probe the internal structure of z ~2-3 galax- 
ies, both in terms of increased flux and source plane res- 
olution. Procuring a larger sample of highly magnified 
z >2 galaxies is thus a highly profitable exercise for fur- 
thering our understanding of galaxy assembly and metal 
enrichment in the Universe. 
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